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On the lattice vibrations in cubic ionic crystals 


By Sric¢ 0. Lunpevist 


SUMMARY 


This paper contains some further considerations connected with a previous paper in which 
the long polar oscillations of a cubic ionic lattice were discussed on the basis of a potential energy, 
which is a first approximation of the expression for the cohesive energy in the Heitler-London 


_ treatment of ionic crystals. A somewhat more accurate estimate of the three-body potential is 


- considered here, giving certain correction terms to the formula for the effective charge derived 


ey 


earlier. The limitations of the Heitler-London method in this connexion are discussed and a 
Semi-empirical approach is suggested including also the polarisation. 


Introduction — 


Most earlier treatments of the vibrations of an ionic lattice have been based on 
the classical theory by Born, according to which the ions are regarded as rigid, non- 


, polarizable structures interacting with central forces. It is now well known, however, 


that the assumptions inherent in the Born model are in general not fulfilled. For 
example, the Cauchy relations are in general not valid in ionic crystals even at zero 
temperature; furthermore it has become clear that the electronic polarisation plays 
an important role; and finally the study of the long polar oscillations has led to 
the introduction of an effective charge, which is known to be about 20 per cent smaller 
than the actual charge for the alkali halides. None of these features can be explained 
on the basis of the original Born model. This paper will present some considerations 
closely connected with a previous paper ([5], hereinafter referred to as A), and it 
will here be shown how the effects mentioned above can be taken into account in 
an approximate way. 

In A a discussion of the vibrations was given, starting from a consideration of the 
expression for the binding energy in the Heitler-London approximation. In ionic 
erystals like NaCl and KCl the overlap between neighbouring ions is rather small 
so that the neglect of third and higher order terms in the overlap integrals is believed 
to be a fair approximation. Neglecting such terms an expression for the binding 
energy is obtained, which besides ordinary two-body interaction terms also contains 
terms representing the interaction between three ions. This is the approximation 
introduced by Landshoff [4] and which was later discussed in much detail by Lowdin 
[7], who first pointed out the occurrence of terms corresponding to many-body 
forces. Léwdin also calculated the elastic constants for several alkali halides taking 
such terms into account. The results show that the many-body forces provide a 
possible explanation of the deviation from the Cauchy relations in ionic crystals, 
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although a discrepancy seems to exist between the observed and calculated value 


of the elastic constants. A different approach is due to Herpin [3], who investigated 
the polarizability forces and found that the quadrupole terms give rise to deviations 
from the Cauchy relations. The calculated deviations are, however, not in agreement 
with those found experimentally at low temperatures, even the sign is incorrect and 


thus it seems as these forces cannot alone be responsible for the deviations at low 
temperatures. The question about the deviations from the Cauchy relations seems not 


settled yet, but is intimately connected with the problem of finding a more accurate 


solution of the quantum mechanical many-body problem for a crystal. For the sake — : 


of simplicity we shall not consider the quadrupole forces investigated by Herpin 


but regard the deviations as due to the:many-body terms occurring in the Heitler- — 


London approximation. 

In a discussion of the long polar oscillations of ionic crystals it was suggested by 
Lyddane, Sachs and Teller [6] that the actual ionic charge should be replaced by 
an effective charge. It was later pointed out by Szigeti [9] that the deviation from the 
actual charge might be due to the mutual distortion of neighbouring ions owing to 
their overlap. In another paper [10] he also derived a macroscopic relation from which 
the effective charge can be calculated from empirical data. In A it was shown that, 
using a very simple approximate expression for the three-body potential appearing 
in the Heitler-London approximation, an effective charge was introduced, which 
could be related to the elastic constants c,, and ¢,,. In this paper a more accurate 
estimate of the three-body potential will be used, giving rise to certain correction 
terms to the formula derived in A. The most obvious way to obtain the effective 
charge in this approximation, however, is to calculate the distortion moment from 
the formula for the charge density." 

It is well known that the electronic polarisation has to be included in order to. 
obtain the correct limiting frequencies for the long polar oscillations and it seems. 
certain that the polarisation will also influence the vibration modes of intermediate: 
and short wavelengths. Unfortunately a satisfactory treatment of the polarisation 
does not exist, but one always assumes that the polarisation of the ions can be treated 
as if the field were homogeneous, having the same value as at the centre of the ion. 
Due to the strong inhomogeneity of the field around a lattice point it must be ques- 
tioned whether such a treatment is at all permissible. For lack of a more satisfactory 
theory it will be used here as a device to obtain the correct limiting frequencies and 
one might at most hope that it might give the correct trend also for waves of inter-. 
mediate and short wavelengths. 


The potential energy 


The discussion in this paper will be restricted to diatomic cubic ionic crystals of 


the NaCl type with ions having complete electron shells, and we shall assume (i) 
that only nearest neighbours overlap and (ii) that the overlap is so small that terms of 
third and higher orders in the overlap integrals can be neglected. For a more detailed 
discussion we refer to A and to the references given there. As in A we divide the poten- 
tial energy of the lattice in the following three parts: 


A. The.classical electrostatic energy of the lattice, which is obtained by considering: 


the ions as point charges 


1 The writer is indebted to Prof. I. Waller for this remark, 
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j °F Efron tol : 
Piero &g denotes the valency of the ion at the lattice point G and e is the-elementary 
4 F Be fe depending only on the distance between two neighbouring nuclei 

f ®, =e > > V’ (\t¢—Te|). (2) 


__ ©. Terms which correspond to a three-body potential. After some elementary 
_ transformations given in A these terms can be written 


Oo=225 > ana £¢ Suy (v@u)e,— 
ite ¥3 
=e > » DEP ee Suv Syu (u Gu). (3) 
B v & 


~Here yu denotes the orbital and spin quantum numbers of an electron associated 
with an ion at a lattice point G,. S,, are the overlap integrals 


Sue = J ph (ae) wy (20) dr — Sup. 


f dz implies integration over the space coordinates and summation over the spin 


coordinates and x means the space and spin coordinates of an electron (x =r, ¢). The 
symbol (yGu) means the integral 


1 
(» Gu) = fe (x) ir opale 1 ke 


(vG 4g, 18 obtained by splitting (v Gz) into two parts by dividing space in two nonover- 
lapping regions containing the lattice points G,, and G, respectively so that 


(vy Gua, + (vGu)e,= (vpn). 


A simple approximate expression for the integrals (y Gu), and (uG) was obtained 


in A by assuming that the integrands are vanishingly small except in a small region 
around G,, so that the denominator |r —r,| can be approximated by |r¢,—Te|. A 


better approximation can be obtained by making use of the expansion of |r —rg|~* 
in spherical harmonics around G,. The expansion contains both positive and negative 
powers of | Te,- Te |, but terms with positive powers will be important only if the points 


G and G,, are so close together that the integrand is appreciable outside a sphere with 
radius | r¢,—I¢|). For crystals with small overlap the wave functions are well localized 
Ub 


so these terms will be fairly small, and we shall for the sake of simplicity omit them 
in our discussion. Taking the first two terms in the expansion into account we obtain 


the following expression 


Oo=P XD O° fa (\re—Terl) 
E ag 


+e >> 2 hal|te-Te'|) 
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The potential energy for small displacements 


As the potential energy (5) is expressed in terms of the relative coordinates of the — 
nuclei it seems most natural to expand it according to powers of the relative displace- 
ments of the nuclei. The treatment of the normal vibrations of a lattice, when a general 
potential is expanded in powers of the relative displacements of the nuclei, has been — 
given by Waller [11] and we shall first review those of his formulae, which will be — 
used in the following. In the expansion of ® the second order term ®, has the form 


O,=3 2 mer y2 Pike ks wie whe : (6) 


Wl" Kk’k’ ap ap 


As usual the index / refers to a lattice cell, & denotes the particles in the cell and ; 
U ik j are the cartesian components of the relative displacements of the nuclei (1, k) and — 
UB). 

‘For crystals where each ion is a centre of symmetry, Waller has found the following 


expression for the elastic constants ; 


1 ree ; 
0 SS LU Ci phy ae \@ 
aypo— 3 be Pree Lick key (7) 
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where A denotes the volume of the elementary cell. It will be of use to have also the 
expansion in terms of the absolute displacements: 
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uantity g¢(req) defined by 
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: a Ae [v4 
c" being the components of r¢¢, we can write (10) as 


~ 


Bes Ga A Ss SO Ss, satee) al é ). (11) 
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om (10) or ay it follows that ®’ vanishes for an undisturbed lattice of the 
"NaCl type. 


The contribution to the coefficients bee” which we shall denote by ¢' ie 


can easily be found be expanding (11). ont we find 


T aye a , ” m” , 
Hef = —e fou (GE TH Aaa + ev [Gk Phe}, (12) 


where we have introduced the notations 


oe ; yf ete tue) er 1 
, Ques leat cuit ane T= | 13 
p be alk, Jake Me Veal only rie No oF 
q Fa ae 

with Gi=0 and Ti =0, 

5 ap ap 

and where (Ge Titer lap = 2 Gc Tie 

a ay yp 


: 


The notation { }, indicates that the derivatives are evaluated at the equilibrium 
SD ictiourdtion: 
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where we have introduced the notation 
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Here ay is the Madelung constant, and a, denotes the interionic distance in equilib- 
; rium. 
By means of (9) we can now find the coefficients Bik. Certain terms drop out 
- , 


because of the symmetry and we get 


Div = —e[e, {ex-+ 12 fx (a9)} fl Oe Le 2Vie Wee ; (16) 
ap ap ap ap ap : 
where ! 
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The effective charge 


In A a formula for the effective charge was derived from the approximate expres- 
sion for the three-body potential considered there. The inclusion of the additional 
term ©’ implies further contributions to the effective charge, which will be derived 
from (16). First, however, we shall derive the effective charge from the definition 
introduced by Szigeti [10] (cf. also Fréhlich [2]) by means of a certain dipole mo- 
ment per unit volume for a small displacement u between the positive and negative 
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; , and I Dein das -eN g, Ng being the number of electrons ~ 
—N, =ég. The contribution from this part of the charge 


y to he ene is immediately found to be a 
= Ba; Pie geen 
The overlap density is hs 
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if we neglect quantities of third and higher orders in the overlap integrals.’ 
- The contribution to the polarisation 


Sexe oe 


Za € 
& P= 75 2 Pw (rly), 
oe a 
V being the volume considered, can be simplified by writing r= Tc, +e and 
- noting that 


3 >, Pur (1 + S)on = Sup = 0 
: “A v 

“and that (ulelm)= 

~ ‘Thus P.= 77 2 Sw (lel): 


| aS (v|p|u) is directed along the vector re,¢,, so writing 


4 eye 
ae it 
v|e| 1) = Mw = 
: Gy Gy, 
. 1 The product wr (x) py (x), w+ is considered as a quantity of the order S)p. 
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A Weer 
denotes a ed ion. Writing Siiye Mm, “ave - ra ha ‘Mz (a) are t 
ments associated with positive and negative ions respectively, we get — 


P,= —A 


I 


Introducing polar coordinates with the origin at the centre of the positive ion and 
the polar axes directed towards the negative ion we get 


M (Ap) a he Pie {2 Sup (v |r cos 6 | LM) — a Spe Soups : (18) 
Lu y p 
Adding P, and P, we find the effective charge 
et = [ BP ac) (a) 
XQ 


or writing e*=e*e 


\ 


* 4m (a) 


ge ee Be (a) (19). 


0 


A relation of this form has recently been derived by Born and Huang [1] from the 
assumption of a local distortion moment acting between neighbouring ions. 
Another formula, which, although formally different, agrees with (19) up to second 


order in the overlap, can be derived from the coefficients Pix (16) in the following ; 


a 
way. We consider two ions (I, k) and U, ) at large distance from each other and 
neglect terms, which tend to zero as (r{i-)° or faster. Then the first term in (16), 


representing a Coulomb interaction, will be retained, View represents the short range — 


forces and gives no contribution. For large rj} it was aioe in A that 
Okie = 4 63, 5 fe (@,) Tie. 
ap ap 
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ae obtain the expression for Wie for large ri. we introduce the explicit ex- 
pression for the quantities G'!,, and get 
: ap 


, 8 hye iin h - , , , 
Wii = ex poem te) +4 (he — ii] i 
aB A ap 


~ Here k’ refers to an ion of the opposite kind as k’. 
For large ri, we get for the coefficients Pin: 


4 re 


8 


ere [ee +12 fir + 4g fie + Bie — Ieee] + — (Ieee — in| Tt, 
0 


aB 


_ which represents a Coulomb interaction with the actual charge replaced by an effec- 
tive charge e* = e*e, where 
8 


SoS ik 
a 


, , 7 3 
Ek lew 12 fee + 4 (dq fe the —hy) + — (he — in | : (20) 


The unsymmetrical appearance of (20) results from the fact that we have neglected 

terms of third and higher orders in the overlap in the formula for the potential energy. 

By calculating the total electrostatic energy of the lattice from the formula for the 

- charge density (17) and deriving the contribution to ¢i%, we should instead have 
aB 


obtained 


, , t 4 
eet 6 fr +2 (a fu + he — hi) + (re — hu) 
0 


* 
Bess 


(21) 


which agrees with (20) up to second order in the overlap integrals (note that the 
quantities f and h are of the second order in S,,,). Although (21) differs in form from 
(19), the two formulae can be shown to be identical by transforming the integrals in 
(21) to a coordinate system with the origin at centre of the positive ion and the polar 
axes directed towards the negative ion. 

The formula given here can be compared numerically with that given by Szigeti 
[9], which is a macroscopic relation and thus independent of any specific model for 
the interaction between the ions. We have computed ¢* for NaCl according to (19) 
from the wave functions for free ions using tables for the radial integrals occurring, 
which were kindly placed at our disposal by Léwdin, and found the value e* = 0.92, 
which is the same value obtained by the approximate formula given in A. The value 
obtained by Szigeti, however, is considerably smaller, he found e* = 0.74 (at room 
temperature). This large discrepancy demonstrates clearly the poorness of the Heitler- 
London approximation. Although the lattice energies and interionic distances agree 
rather well with the empirical values, this is not the case with other properties, which 
depend more sensitively on the wave functions used. For that reason a calculation 
of the lattice vibrations cannot be successfully done of the basis on the energy terms 
calculated in the Heitler-London approximation. It seems therefore reasonable to try 
to determine the parameters appearing in such a way that the most important 
quantities like the interionic distance, the elastic constants and the effective charge 
will agree with experiment. It would of course be most desirable to consider interac- 
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ance and by using the equilibrium condition. The remaining pe 


troduce some reasonable assumption as regards the overlap 


in order to assign definite 
values for these parameters. 4 
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Approximate treatment of the polarisation 


It has been shown that the deviations from the Cauchy relations and the introduc- 


tion of an effective charge are natural consequences of a quantum mechanical — 


be chosen so as to give the correct effective charge, but in addition one has to in- 


treatment according to the Heitler-London approximation, even though the numerical — 


agreement is not very good. Nevertheless the treatment given here does not provide . 


a satisfactory starting point for calculating the lattice vibrations, the main reason 
being that the Heitler-London approximation does not account for the electronic 
polarisation. From the study of the long polar vibrations it is well known that the 
neglect of the polarisation leads to incorrect limiting frequencies and it seems probable 
that also vibration modes of short and intermediate wavelength will be strongly 
influenced. Unfortunately a satisfactory treatment of the polarisation of ions in a 
lattice has not yet been given, but one always assumes that the polarisation can be 


= th tlt 


calculated as if the field were homogeneous. Due to the finite extension of the ions | 


and the strong inhomogeneity of the field this assumption seems far from justified, 


but nevertheless we shall make use of this assumption here and derive the contribu- — 
tion to the coefficients i. A treatment for a general lattice has been given by Born’ 


a 
and Huang [1], who conaiier waves and apply Ewalds transformation to the dipole 
field. Here we shall confine ourselves to a lattice of the NaCl type and by first consider- 
ing a large but finite crystal instead of an infinite one the formal treatment can be 
made simpler. 
The field at (J, ) due to a system of dipoles pj, associated with the other particles 
(U, k’) is 
Ey= > D> The Pie. (22) 
og Vk’ B ap B 


We shall now calculate the exciting field, that is the field at the instantaneous 
position of (/, k) due to all the other particles. This field is the sum of three parts: 
(i) The Coulomb field corresponding to the displacement of the nuclei I’, k’), which 
is obtained by putting ph = ex uy. 
(ii) The field due to the induced dipoles: pi. = wie 


(iii) The field change owing to the displacement of the particle (I, k) itself, which | 


can be calculated as if the particle is kept fixed and that all other particles are 
subjected to a displacement — uj, thus Pe = en-us. 


Denoting the exciting field by Ej, we then obtain 
Br= >, > The (cheat. + ule). 
ie Se wes Rica ial er 
We introduce the notation 
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rs E and E°, we get 
E-E)+Tak. (25) 


~ 
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If det (1 — Ta)+0, which we shall assume, we get the formal solution 


ag = Balbaet FE. (26) 
mM Yow the force on the aceeas (i, k) uN to the field is simply e% Ei,. This includes the 
E classical coulomb field as well as the coulomb part of the many-body interactions, 
_ which we have already considered. The additional force due to the polarisation is 
found by subtracting the field E°, thus 


Ki= et {T QE}, (21) 
a a 
where we have introduced the symbol 
a a 
- gt ' 
; Q= (5-1 27" 


Biriderting the eps for the coulomb field E° we can easily find from (27) the 
~ coefficients p- @, associated with the polarisation forces 
ap 


Boke = ek Ch {T Q Thee. (28) 
ap ap 


Neglecting surface effects these quantities will not depend on / and I’, but only on 


the difference /-I’. 
The polarisation forces can alternatively be derived from the following expression 


for the potential energy 
Oo? = — E°QE’. (29) 
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which can easily be reduced to the standard form (6), giving the coefficients 


; “rie = — eh hen {Dice [Q Tee tap- (30) 
a 


The coefficients $7%.- can be calculated according to (9) and noting that sums 


: B ; 
of the type > ex Te, vanish, we immediately get (28). 
Vk’ ap 


Inserting (30) in the formula for the elastic constants (7) we find that the polarisa- 
tion forces do not contribute to the elastic constants. 


The discussion given here has been based on the formal solution (26), but the coeffi- _ 


cients Pix can hardly be computed according to (28) because of the difficulties in 
a 


performing the matrix inversions and multiplications. In order to investigate the 
lattice vibrations, however, we only need the associated Fourier matrix. This is 
fortunately rather easy to calculate, for by making use of the faltung property the 
calculation will be reduced to operations on 6 x 6 matrices, which can be handled 
by standard methods. 

In A the polarisation was included for the long optical vibrations by means of 
certain macroscopic relations involving the dielectric constants for static and high- 
frequency fields, and formulae for the transverse and longitudinal frequencies were 
derived, which were generalisations of those given independently by Odelevski [8] 
and Szigeti [10] for the case of central forces. These formulae can now be obtained 
from (16) and (28) by neglecting the non-coulomb part of the many-body terms in 
(16) and by making use of the relations between the elastic and dielectric constants 
and the parameters in (16) and (28). A corresponding derivation has been given 
by Born and Huang [1] in the case of central forces, so the details will not be given 
here. 


Comments 


In this paper a set of formulae have been derived from which the lattice vibrations 
of ionic lattices of the NaCl type can be calculated, allowing for the effects of the 
deviations from the Cauchy relations, the effective charge and the electronic polarisa- 
tion. An accurate treatment of this problem would require an accurate solution 
of the quantum mechanical many-body problem for a crystal and therefore a semi- 
empirical approach has been suggested here, based mainly upon the Heitler-London 
approximation. The Heitler-London approximation leads without extra assumptions 
to deviations from the Cauchy relations of the correct sign and order of magnitude 
and to the introduction of an effective charge, although the form of the potential 
considered here is certainly too restrictive. From several points of view, however, 
the approach suggested here is not very satisfactory. The terms in the potential 
energy, which are associated with the deviations from the Cauchy relations and the 
effective charge, are derived from a very approximate quantum mechanical treatment 
and the polarisation has been treated on a classical model. It seems not at all obvious 
that the potential energy derived from an accurate quantum mechanical treatment 
of the crystal could to a fair degree of approximation be reduced to the simple form 
discussed here. Therefore, the treatment given here ought to be considered more as 
a phenomenological theory, the aim of which is only to take into account certain 
features not considered in the classical Born model. 
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